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218. 


A THIRD MEMOIR ON THE PROBLEM OF DISTURBED: ELLIPTIC 
MOTION. 


[From the Memoirs of the Royal Astronomical Society, vol xxxi, pp. 43—56. 
Read January 10, 1863.] 


THE object of the present Memoir is to obtain the differential equations for deter- 
mining 
r, the radius vector, 


v, the longitude, 


y, the latitude, 


of the disturbed body, when the last two coordinates are measured in respect to an 
arbitrarily varying plane (which however, to fix the ideas, is called the variable ecliptic) 
and the departure point or origin of longitudes therein. This is very readily effected 
by means of an expression for the Vis Viva function given in my “Supplementary 
Memoir on the Problem of Disturbed Elliptic Motion," Mem. Roy. Ast. Soc, t. XXVIII pp. 
217—234 (1859), [215]. Neglecting the squares of the variations of the variable ecliptic, 
and also the products of the variations by sin y, or zd then (as might be expected) 
it is found that the equations for r and v are the same as for a fixed ecliptic, and 
the equation for y is found in a simple form, which is ultimately reduced to coincide 
with that obtained for the lunar theory by Laplace in the seventh book of the 
Mécanique Céleste, and which is used by him to show that the effect of the variation 
of the ecliptic on the latitude of the Moon (as measured from the variable ecliptic) 
is insensible. And it is shown conversely how the approximate formula of the Memoir 
may be obtained by a process similar to that made use of in the Mécanique Céleste. 
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L 


The position in respect to a fixed plane of reference and origin of longitudes 
therein, of the variable ecliptic and of the departure point or origin of longitudes 
therein, are determined by 

0', the longitude of node, 


o’, the departure of node, 

¢’, the inclination, 
where, by the definition of a departure point, 

do’ — cos $' dO’ = 0, 

and then, in respect to the variable ecliptic and departure point or origin of lougi- 
tudes therein, the position of the disturbed body is determined by 

r, the radius vector, 

v, the longitude, 

y, the latitude; 


and this being so, then (Supplementary Memoir, pp. 220, 227) the expression for the 
Vis Viva function is, 


T=} [e +r (Œ +R}, 
where : 
Q= S + [eos (v — e^) sin $' . & — sin (v — o^) ġ'], 


R= cos y . 9 — sin y [sin (v — o^) sin œ’. & + cos (v — e^) $7, 


the superscript dots being used to denote differentiation with respect to the time. 
The last-mentioned expressions may for shortness be denoted by 


Qiz —y+t+A, 
R=cosy.s—Bsin y. 


The equations of motion are of course, 


d dT dT dV 
dt dr dr dr’ 
d dT dT dV 
dt dv dv adv’ 
d dT JU a. 


d dj^ dy" dy’ 
2693 
where V= ne +Q, if Q is the disturbing function, taken with Lagrange’s sig? 


(Q=—R, if E is the disturbing function of the Mécanique Céleste). 
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To reduce these, we have in the first place 


dT 2 aT 
dr LÀ , gr 77 (Q1, 


dT dT dQ : 
qj ^" Rosy, Fh =n (Qe +Ro)= r (—QB-—RAsiny); 


oF NO ‘ nek =r R(—siny.t—Becosy). 


The equations are thus reduced to 


dr dQ 
dt? -r(Qe Rer pem 


Í (22008 y . B) +r (QB+RAsiny) = 5- 


d "M, dQ 
5 C79 ) trEQGhy; ra 
and then substituting for Q, R their values, viz. 


d 
Q= -t4 


R= cosy — B sin y, 


we find 


d?r dv (dyV| ., n'a? dQ ; 

dB r Joos v (Fe) +(e - ree 

d dv i" 

di (recos y F) sae +8, 
y n,(*y a 

d (rH) +r cos y sin y (T) =W 


d : dv , 
9 -r(-24 4Y _ 9B sin y cos y Gp +A? + B sin’ y) 


+ 6, 


where 


= (r Bsin y cosy)  * (5 - 4 sin y cosy 5 — AB cos y ), 


6 =F A) +r (- (cos? y — sin” y) Ba at + B'sinyeosy), 


in which de dej 
A = cos (v — c") sin $^ 5, + sin (v — e") 24 , 


i 20. n ad’ 
B — sin (v — o^) sing W + 008 (0 — 0) 5 


6’, c', d/ being given functions of t such that do’ —cos ¢' dd’ = 0. 


The foregoing equations of motion are rigorously accurate. 


64—2 
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II. 
Neglecting the terms which involve A?, AB, B*, we have 


A = r(- 24 Y - 2B sin y cosy $9), 


In x dy F dv 
= G(r? Bsiny cosy) +r (B — Asin yeos y 5) ; 


d 
=p 4) — r? (cos? y — sin? Bt, 


and if we then neglect also the products of A and B by y and A we have 


where it may be noticed that, in order to obtain this last value of ©, the only 
neglected term is a term containing B sin? y. 
Now, attending to the values of A and B, we have 
dA dA dv dA 
dt dv dt dt 
dv | dA 
di* dé’ 


where here and in the sequel 3 denotes differentiation in regard to ¢ in so far only 


=—B— 


as it enters through the quantities o’, 0’, ¢’, which determine the position of the 
variable ecliptic. 


Hence * 
dr dv dA v 
¢= 2Ar a+ e(- Bo +S) BS 
dr d A dv 
= 2Ar 7 tU dr di — 2B; 


and, as above, W= 0, $ =0. 
Let r, v y be the values obtained on the supposition that ($—0, and 
r 4-Ór, v+ dv, y + by, 
the accurate values; the first and second equations show that, neglecting the products 
of y and es into dy and d we have 5r=0, 8v—0; so that the values of r and " 


are not affected by the variation of the ecliptic. And then, substituting in the third 
equation y+ôy in the place of y, and for 


cos (y+ dy) sin (y+ êy), = cos y sin y + (cos? y — sin? y) dy, 
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writing cos y sin y + dy, we have 

S0 es e) (cos ysiny +8) =F +8 ST tS 

or, since the terms independent of dy and © must destroy each other, this is 
See (FY Sy = um at 

or, as this may be written, 


d?ày dr dày dv dQ, dr d'A dv 
10 art E +e) yc rgt E B 


dt dt dt ' 
that is, 


addy 12 dr dày m) Bs 1,.dO ldr dA dv 
df rdi di ie 


pee ior a a ge 28 


or, what is the same thing, 


dS 1 dQ P n dv\? 

"E te + hm) by = 5 BG + ee mene — (2) bby 
2 dr dày l dr dv 
rdi di $24 Gt 2B (n x3 
EY M. 


, 


where as usual m, =“, is the ratio of the mean motion of the Sun to that of the 


Moon; the term 3 mn? ôy having been added on each side of the equation in order 


dQ, 


to destroy on the right-hand side the corresponding term arising from 18 dp In 


' ; Q 
fact, to find the approximate expression for T , we have 
^y 
= T. cos? H — 4), 


where H is the angular distance of the Sun and Moon; that is, cos H = cos y cos (v — v) ; 
here 


mr? 
Q = = {3 cos? y cos? (v — v^) — $}, 


OF T, 3 sin y cos y cost (v — v) 


mr 


don "(3 +3 cos (2v — 2v’)) sin y cos y ; 
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or, neglecting the periodic quantity cos (2v—2v), and writing y for siny cosy; also 
putting as usual m = n’a” =m na", and 1 =a’, we have 


dQ ^ 

dor NT 
and thence ue 

óc o5 dx 
id l 0 


Substituting this value of = TES , and putting also r=a, To 0, o 0, the 


equation for dy becomes 


d*à dA 

ae cU + 3m?) dy = — 2Bn + dr 
III. 

To deduce the formula, seventh book of the Mécanique Céleste, I proceed as follows: 

Putting 

dóü 1 dec 

dt cos dt’ 


we have 


A cos ¢’ = do. sin d Wc. cos $' sin (v — a’), 


B cos ¢’ = ae sin $' sin (v — o’)+ i cos $’ cos (v — a’), 


which may be written 


A cos $* - — sin v © (sing cos o^) -- cosu $ , (in g sin c^), 


d STU 
B cos$'— cos "d ,(sin Q' cos o^) 4- sin vd © (sin ¢ sino’). 


Laplace in effect assumes that the variations of the ecliptic are given in the form 
sin $' sin o’ = — Xk sin (int + €), 
sin g’coso’= Èk cos (int + e), 


(it+e is there written for the argument, n being assumed =1) where 2, k, e are 
absolute constants, the quantities 4 being all very small in comparison with m? Sub 
stituting these values, and putting cos ¢’ equal to unity, we have 


A — — ik cos (v + int + €), 
B = — X4k sin (v + int + e); 
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and thence also 
—= n> řk sin (v+ int +e), 
and 


nB — 1 — — nx (2i +ë) ksin (v+ int +6), 


so that the equation for dy becomes 


d?ày 
dt? 


+n? (14-8 m?) dy -- nX (2i 4-77) k sin (v 4- int + e) 20; 


or, taking as the independent variable v(= nt) in the place of t, this is 


d*ày 


2 E + § mt) 8y +23 (2i P) ksin(v+iv+e)=0; 


which is, in fact, Laplace’s equation, n being retained instead of being put equal to 
unity, and dy being the part which depends on the variation of the ecliptic, of his s,. 


EV. 
Conversely the equation 


TM nt (1 +3 m?) dy = E E 


may be obtained by a process similar to Laplace’s. Assuming that the Moon and 
Sun are each of them referred to a fixed plane of reference and origin of longitudes 
therein, by the coordinates 

u, the reciprocal of the reduced radius vector, 

v, the longitude, 

s, the tangent of the latitude, 


for the Moon, and by the corresponding coordinates w’, v', s’ for the Sun, then we have 


di dO . dQ 
ds dedo ^^ BOG d 
Rec TIE n 
nu (1+ 2. fag, d") 


Here, as before, 


a=% qe H- 9, 


or, as 1t is now to be written, 


m ' (14-8) w^? (s (v — v) =) " i ] 
"sse ÜWrrevrrss CJ) 
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or, since the second term 
mu’ plte 
(1+ gay i 


gives, as is immediately seen, no term in 


ds dQ dQ dQ 
do dv du (t9: 


we may, in calculating this quantity, write 


s m (1lo-5)w* , (= (v — 7) 4- =A 
(+s) w P Wi4+sV14 sg 
4 2 
23 mee p (cos (v — v) 4- ss’ . 
or, neglecting s^, 
Q— mu’ (= (v— v) + S . 
u 


Hence, putting for a moment 


cos (v — v^) + ss' 
6 z————————, 
u 
we have 


dQ dQ jc ds dO de 


ds dQ 


de 
dv dv ate it 


where the factor in { } is 


= 7 {- SF sin(o— 0) +8 (cos (wv) +e) - a +s) 


1 "AP Iur Rand 
=> [es (v =v’) — sinu- , 
and the above-mentioned quantity is 


3m^u^? 


^ - V MM ge M 
-—— (gos o — enr fs cos w- 1) -q sin n f) s}; 


2 


or, what is the same thing, 


3m'/u* 


cos (v — v’) reos (v — v) F sin w=}. 


Considering now the Sun as moving in the variable ecliptic, its latitude is 
= sin $' sin (v' — o^); that is, we have 


s' — sin $' sin (v — e); 
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and if the Moon were in the variable ecliptic, its latitude would be sin$'sin(v— o"); 
that is, the latitude, measured from the variable ecliptic, is —s— sin ¢’sin (v—o’); or, 
putting 

s,, the Moon’s latitude, measured from the variable ecliptic, 


we have 
s=sin $' sin (v — e^) 4- s,. 


Hence, disregarding the variations of $', o’, we have 


d. 


dE. mu S. 
dy" né cos (vy — 7) + 57 : 


and substituting these values of s, E , and s', we find 


sos (o —v) —$ sin (v— v) - s 


ds, 


dv 


= s$,c08(v — v') — 7^ sin (v — v^) 
+ sin $ [sin (v — a^) cos (v — v’) — cos (v — c^) sin (v — v) — sin (v — 7) 


=  8,C08 (v — v) - sin (v - 9); 


or 


eur We (v — v) fs, cos (v — 9) — e sin (v — v) ; 


which, neglecting the periodic terms, is 


m'u 
a $ y? S, > 
and then 
ds dQ, dO ,dQ 
dv du T du T F9) qon g mu” Py mU” ; 
naa iwi ien 
a h? J udv 
: eh 1 
which, putting m =n” a^? = mn? a^, u "a P qam D becomes 
- imis; 
so that the differential equation is reduced to 
2 
TS bet} ms, = 0. 
C. III. 65 
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But 
s = sing’sin(v—o’)+s, 
os = sng’ cos (v — o) + © sin $ sin(v—o) +, 
d?s "I A up ee Lone 8 d's, 
ae sin $' sin (v o’) + sin $ cos (v — 0) + 7; 2: sin $' sin (v — o^) + 5 d: 
where : denotes differentiation in respect to v, in so far only as it enters through 


$' and o’ (these are functions of f, which is -1 »). Hence 


d?s ds, d d d 
dp *? “ah * sin $' cos (v — e) t5, 3,9 9 sin (v — e^); 
but 
" PR ee AA TAN. E ; E peters Ving 
gr d sin (v — 9) = ~ (sin v 4; sin $ cos a! — cos v 7, sin d'sin a’) 
1 ; 
=- z 4 cos ġ 
n 
and thence 


LAA 1dA 1 1 d'A 


d d' 
dv dv sini) en(epr)—- ee de n mdi 


and similarly 


d o " 1 di s ox TAE UENIRE REN 
g 9in $ cos (v7 0) = T feos v 5, (in d cosa’) +sin v 5 sin ġ sin o'} 


Hence 
d?s ds, .2 l1dA. 
ip tt ga MOT GS ero. 


or, putting v = nt, 
d?s 1 d’s 2 1d A. 
n 


de mde tet, 3m a E 
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whence, substituting in the equation 


2 
TS sd més, =0, 
we have 
Qo à CA... 
qon (1 +m)s, +2nB--7 =0; 


or putting s, — y + dy, then 


d? 
atm rim)y- (0 


which gives, in the approximation which is being considered, the principal term of the 
latitude, and then 


d'ày 


apt n? (Lt 10d) Sane ue 0, 


dt 


which is the approximate equation previously obtained by the method of the present 
memoir. 


for) 

Qı 
| 

[59] 
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